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Abstract 

In this paper we perform an asymptotic analysis for two different vanishing 
viscosity coefficients occurring in a phase held system of Cahn-Hilliard type that was 
recently introduced in order to approximate a tumor growth model. In particular, 
we extend some recent results obtained in [5], letting the two positive viscosity 
parameters tend to zero independently from each other and weakening the conditions 
on the initial data in such a way as to maintain the nonlinearities of the PDE system 
as general as possible. Finally, under proper growth conditions on the interaction 
potential, we prove an error estimate leading also to the uniqueness result for the 
limit system. 
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1 Introduction 


In this paper we study the system of partial differential equations 

adtii + d t tp - Ap = p(<p)(a - 7 / 1 ) (1.1) 

p = /3d t tp - A ip + F'(ip) (1.2) 

d t a - Act = -p((p)(a - 7 / 1 ), (1.3) 

together with the boundary and initial conditions 

d u p = d u tp = d u a = 0 (1.4) 

/i(0)=/i 0 , <p(0) = <p 0 and cr(0) = cr 0 . (1.5) 


Each of the partial differential equations (j 1.1 [i — (1 1.3 j) is meant to hold in a three-dimensional 
bounded domain f 2 , endowed with a smooth boundary T, and for every positive time, and 
d v in (II.4p stands for the outward normal derivative on T. Moreover, a and (3 are non¬ 
negative parameters, strictly positive in principle, while 7 is a strictly positive constant. 
Furthermore, p is a nonnegative function, and F is a nonnegative potential. Finally, /i 0 , 
<Po and (To are given initial data defined in 

The physical context of this paper is that of tumor growth dynamics. This topic has 
in recent years become of big interest in applied mathematics, especially after continuum 
models were developed (cf., e.g., [71H8] h The fact that multiple constituents interact with 
each other made it necessary to consider diffuse interface models based on continuum 
mixture theory (cf., e.g., [4Mlll2M15M20M26| h These models consist of a Cahn-Hilliard type 
equation (in general with transport) containing reaction terms that depend on the nutrient 
concentration (e.g. oxygen) and in turn obey an advection-reaction-diffusion equation. 
Even though numerical simulations of these models have already been carried out in 
several papers (cf., e.g., [7J Chapter 8 ] and references therein), the rigorous mathematical 
analysis of the resulting PDEs systems is still very poor. To our knowledge, the Erst 
results are related to the so-called Cahn-Hilliard-Hele-Shaw system (cf., e.g., [TUI [25]) 
in which the nutrient is neglected, while two very recent contributions [5] and [133 (cf. 
also [16], where formal studies on the corresponding sharp interface limits are performed) 
deal with a model recently proposed in m (or approximations thereof, see also |23), 
where the velocities are set to zero and the state variables are the tumor fraction <p and 
the nutrient-rich concentration < 7 . We can set <p ~ 1 in the tumorous phase and <p ~ — 1 
in the healthy cell phase, while a typically satisfies a ~ 1 in a nutrient-rich extracellular 
water phase and a ~ 0 in a nutrient-poor extracellular water phase. Moreover, the third 
unknown p is the related chemical potential, specified by (j 1.2 [) as in the case of the 
viscous Cahn-Hilliard or Cahn-Hilliard equation, depending on whether (3 > 0 or /? = 0 
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(see mum). In addition, in [2] the PDE system (11.10 - (11.511 was studied for the very 
particular case that a = f3, and the asymptotic analysis as the coefficient a = f3 tends 
to zero was performed, yielding the convergence of subsequences to weak solutions of the 
limit problem; moreover, in ra the existence of weak solutions, as well as uniqueness and 
existence of attractors, was proved directly for the limit system where a = (3 = 0 (cf. also 
the following comments in this Introduction). 

In the case a = 0, the sub-system (11. ID — (11.211 becomes of viscous or pure Cahn-Hilliard 
type, depending on whether /3 > 0 or j3 = 0. On the other hand, in the case a > 0 the 
presence of the term ad t p in (11.111 gives a parabolic structure to equation (11.11) with 
respect to /i. 

We remark that the original model deals with functions F and p that are precisely 
related to each other. Namely, we have 

P(u) = 2p 0 y/F{ u ) if |u| < 1 and p(u) = 0 otherwise, (1.6) 

where p 0 is a positive constant and F(u) is the classical Cahn-Hilliard double-well free 
energy density. However, this relation is useless in many aspects of the mathematical 
study. Moreover, one can allow F to be even a singular potential. 

As mentioned above, [5j just deals with the case a — [3 for the mathematical study, 
although the constants a and (3 have a different meaning. In that paper, the existence of 
a unique solution to the system (11.ID (11.511 was proved under very general conditions on p 
and F, and, in the same framework, the long-time behavior of the solution was discussed. 
In addition, in a more restricted setting for the double-well potential F, [5| investigated 
the asymptotic behavior of the problem as the coefficient a = /3 tends to zero, finding the 
convergence of subsequences to weak solutions of the limit problem. Moreover, under a 
smoothness condition on the initial values, uniqueness for the limit problem was proved 
and, consequently, also the convergence of the entire family. It must be pointed out that 
a uniqueness result was proved in p2J under weaker assumptions. 

In the present paper, we first extend some of the results of [5]. Namely, we let the posi¬ 
tive parameters a and (3 be independent from each other, and we weaken the assumptions 
on the initial data while keeping the potential as general as possible. At the same time, 
we establish a general a priori estimate that is uniform with respect to the parameters 
a and (3. This is the starting point of possible asympotic analyses with respect to these 
parameters. Then, we confine ourselves to a class of regular potentials. In this framework, 
we state a convergence result as both a and /3 tend to zero independently, and we prove 
an error estimate in terms of a and /3 for the difference of the solution to (ll.ip - fll.5j) and 
the one of the limit problem. The case of just one of the parameters tending to zero is 
the subject of a work in progress. 

Let us express our belief that the results of the present paper are general and interesting 
enough so that methods and estimates could be extended to other situations. In particular, 
in case of the trivial choice p = 0 (admitted by our assumption (12.31) ) our system (II.II) - 
(11.31) decouples and (II. II) - (11.21) reduces to a well-known phase field system of Caginalp 
tye which can be seen as a (doubly) viscous approximation of the Cahn-Hilliard system 
obtained at the limit as a and /3 go to zero. To this concern, let us quote the papers 
[81191122023] , where different investigations on this kind of viscous approximations of Cahn- 
Hilliard system are performed, and point out that the results contained in [22] are here 
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generalized and somehow improved. 

Our paper is organized as follows. In the next section, we will state the assumptions 
and our results on the mathematical problem. In Section [3j we will prove the extensions 
mentioned above. The last section is devoted to the asymptotic analysis and the error 
estimate. In the remainder of the paper, we take 7 = 1, without loss of generality. 

2 Statement of the problem and results 

In this section, we make precise assumptions and state our results. As in the Introduction, 
R C M 3 denotes the domain where the evolution takes place and T is its boundary. We 
assume R to be open, bounded, and connected, and T to be smooth. Moreover, the symbol 
d v denotes the outward normal derivative on T. Given a final time T, we set 

Q := Q x (0, T ) and E:=Tx (0, T ). (2-1) 

Moreover, we set for brevity 

V := H\n), H : = L 2 (R), and W : = {v G R 2 (R) : d v v = 0 on T}, (2.2) 

and endow these spaces with their standard norms. For the norm in a generic Banach 
space X (or a power of it), we use the symbol || • \\x with the following exceptions: we 
simply write || • || p and || • ||* if X = L P (Q) or X = L P (Q ) for p G [1, +oo] and X = V*, 
the dual space of V, respectively. Finally, it is understood that H is embedded in V* in 
the usual way, i.e., such that (u,v) = f n u v for every u G H and v G V, where (•, •) 
stands for the duality pairing between V* and V. 

As far as the structure of the system is concerned, we are given two constants a and 


f3 and three functions p, B and n satisfying the conditions listed below 

a,0G(O,l) (2.3) 

p : M —* [0, Too) is bounded and Lipschitz continuous (2.4) 

B : R —* [0, Too] is convex, proper, lower semicontinuous (2.5) 

7 f G C 1 (M) is nonnegative, and 7r := 7r' is Lipschitz continuous. (2.6) 

We also define the potential F : M —y [0, Too] and the graph B in M x R by 

F := B T 7f and B := dB. (2.7) 


We notice that if F is a C 2 function then our assumptions imply that F" is bounded from 
below. We also remark that B is maximal monotone. In the following, we write D{B) 
and D(B) for the effective domains of B and R, respectively, and we use the same symbol 
B for the maximal monotone operators induced on 1 2 spaces. 

Remark 2.1. We notice that, among many others, the most important and typical exam¬ 
ples of potentials fit our assumptions. Namely, we can take as F the classical double-well 
potential and as the logarithmic potential, which are defined by 

F ci (r) := J(r 2 - l) 2 = J((r 2 - 1) + ) 2 + 1((1 - r 2 ) + ) 2 for re* 

Fiagiv) := (1 — r) ln(l — r) + (1 + r) ln( 1 + r) + k(1 — r 2 ) + for |r| < 1, 


( 2 . 8 ) 

(2.9) 
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where the decomposition F = B + 7 ? as in (12.7(i is written explicitly. In (12.9]) , a is a 
positive constant which, depending on its value, does or does not provide a double well, 
and the definition of the logarithmic part of Fi og is extended by continuity to ±1 and by 
+00 outside [—1,1]. Moreover, another possible choice is 

F(r) := J(r) + ((1 — r 2 ) + ) 2 for r G M, (2.10) 

where I is the indicator function of [—1,1], which takes the value 0 in [—1,1] and +00 
elsewhere. For such an irregular potential, the associated subdifferential is multi-valued, 
and the precise statement of problem (11.ID — ('ll,5j) has to introduce a selection £ of B{u). 


As far as the initial data of our problem are concerned, we assume that 

yfa /i 0 , uo e H, tfo G V, and F(ip 0 ) e L 1 ^), (2.11) 

while the regularity properties postulated for the solution are the following: 

H,<T G H\0,T;V*) nL 2 (0,T;R) (2.12) 

<p G H\0,T;H) nL 2 (0,T;IF) (2.13) 

( 6 L 2 (0,T; F), and £ E B(u) a.e. in Q. (2-14) 


We notice that (I2.12j) - (l2.13j) imply that / 1 , a G C°([0,T];R) and ip G C°([0, T]; F). At 
this point, we consider the problem of finding a quadruplet (/x, </?, cr, £) with the above 
regularity in order that (/ 1 , </?, cr, £) and the related function 


R = p((p)(cr - v) 


(2.15) 


satisfy the system 


a{d t fjL,v)+ / d t <pv 


V/i ■ Vu = / Rv 


for every v G V, a.e. in (0, T ) 

H = f3d t p> — Atp + £ + 7r(</?) and £ G B(ip) a.e. in Q 


(•d t a , v) 


Rv 


' V ct ■ Vn = — 

n Jn 

for every v G V, a.e. in (0, T ) 

MO) = ho, <M°) = Fo and cr(0) = a 0 ■ 


(2.16) 

(2.17) 

(2.18) 
(2.19) 


This is a weak formulation of the boundary value problem (11. lj) (1 1.5 ji described in the 
Introduction. The homogeneous Neumann boundary condition for ip is contained in (12.13jl 
(see (12.21) for the definition of IF), while the analogous ones for /1 and cr are meant in a 
generalized sense through the variational equations (12.161) and (j2.18j) . We notice once and 
for all that the addition of (I2.16|) and (I2.18|) yields 


(d t (afi + ip + a),v) + / V(/i + a) ■ Vv = 0 

Jn 

for every v G V, a.e. in (0 ,T). We also set for convenience 

s = Vp(f) M - h) • 


( 2 . 20 ) 


( 2 . 21 ) 


Our first results deal with the well-posedness of the above problem and general a priori 
estimates. Namely, we have: 
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Theorem 2.2. Assume (12.3ft - (12.71) and (12.lift . Then, for every a, /3 6 (0,1), there exists 
a unique quadruplet (ii,<p,cr,£) satisfying (12.12ft - 02.14ft and solving problem (12.15ft - 02.19ft . 

Theorem 2.3. Assume (12.3ft - (12.71) and (12.11ft . Then, for some constant C that depends 
only on R, T and the shapes of n and p, the following is true: for every a, (3 £ (0,1), the 
solution (/i, <^,cx, £) to problem (12.15ft - (12.191) with the regularity specified by (12.12ft - (12.14ft 
satisfies 

tt 1 ^ 2 ||A t IU oo ( 0 ,T;H) + II VT\\l 2 ( 0 ,T;H) 

+ /3 1/2 \\dtT\\L 2 (0,T-,H) + |M|l°°( 0,T;V) + ||-F(</?) |lLi( 0i r;Ll(Q)) 

+ H cr ||iI 1 (0T;E*)nL° o (0,T;//)nL 2 (0,T;I/) + II ‘S'llZ, 2 (0,CT ; J/) + \\R\\l 2 (0,T;H) 

+ || d t (afa + <^)||l 2 (o,t ; v*) 

< C (a 1 ^||/i 0 ||h + 11(^0II v + l|-f 1 ( < ^o)||]/i 2 (Q) + Ikolln) (2.22) 

as well as 

IHI L 2 (0,T;V) + IMU 2 (0,T;WO + || f || L 2 (0,T;/I) 

< C (a 1 ^ ||/io||// + ||(^oII v + II^(^o)llzA*(^) ll <T °H H IMU 2 (o,t-h) + l)- (2.23) 

Thus, a uniform estimate for the left-hand side of (12.221) holds in terms of the norms of 
the initial data related to (12.11ft . while an estimate for the left-hand side of (I2.23|) follows 
whenever a bound for ||/i||L 2 (o,T ; r/) has been proved. 

Remark 2.4. We note that Theorems 12.21 and 12.31 improve the results of [5], since the 
stronger assumption made there, 

ho, To, or 0 e V and F((p 0 ) e 1/(0), (2.24) 

is now replaced by (12.111) (and also since just the case a = (3 is dealt with in j[5j). 

Our next results regard the asymptotic analysis as the coefficients a and f3 tend to 
zero, independently. To this end, we restrict ourselves to a particular class of potentials. 
Namely, we also assume that 

D(B) = R. and \B°(r)\ < C ^R(r) + lj for every r G M, (2.25) 

where B° is the element of B with minimal norm and C is a given positive constant. 
Let us note that, for example, all polynomially growing potentials, as well as exponen¬ 
tial functions, comply with our assumption (12.251) . Let us point out that (12.25ft implies 
(actually, it is equivalent to) the condition 

D(B) = K, |s| < C ^ B{r ) + lj for all r G K, s 6 B{r) (2.26) 

for the same constant C , as checked precisely in the next remark. 

Remark 2.5. In fact, a similar equivalence holds for a more general growth condition 
and in the general setting of Hilbert spaces, as we show at once. If X is a Hilbert space, 
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B : X — >■ [0,+oo) is convex and l.s.c. (thus continuous since it is everywhere defined), 
B := dB and, for every u G X, B°(u ) is the element of B(u) having minimal norm, the 
assumption 

\\B°(u)\\ x < T( B(u )) for every u G X, 
where 'h : [0, +oo) —>■ [0, +oo) is continuous, implies 

IICIIx — f° r every u G X and every ( G B(u). 

Indeed, for arbitrary u G X , ( G B[u) and £ > 0, we have 

(B°(u + £() - C, (« + eC) - u ) > 0, whence ||C|| X < ||-B°(w + eC)|| x . (2.27) 

By applying our assumption to u + e(, we deduce that 

lieu* <*(£(«+<))• 

By taking e —>■ 0 and owing to the continuity of o B, we conclude. 

Now we are ready to state our result on asymptotics. 

Theorem 2.6. Assume (12.311 - 02.7H and (12.251) on the structure and (12. 1 1)1 on the initial 
data. Moreover, let (ha,/3, Pa,/3, &a,/3, £a,p) be the unique solution to problem (12.15p - fl2.19f) 
given by Theorem \2.2l Then, we have that there exists a quadruplet such that 


ha ,/3 —> h weakly in L 2 { 0, T ; V) (2.28) 

Ta,/3 T weakly star in L°°(0,T;R) fl L 2 (0,T]W) (2.29) 

a a ,p —> cr weakly in H l (0,T-,V*) D L 2 (0,T;R) (2.30) 

dt{a ha, p + Ta,p) -t dt<p weakly in L 2 (0, T; V*) (2.31) 

ia ,/3 —> £ weakly in L 2 (0,T;i7) (2.32) 

at least for a subsequence. Moreover, every limiting quadruplet (n,<p,a,£) satisfies 

(d t (p,v) + J n V/i-Vu = f n Bv VvgV, a.e. in (0, T) (2.33) 

h = — Xip + £ + 7r(</?), £ G B(ip) a.e. in Q (2.34) 

(d t a,v) + f Q Vcr ■ Vv = — f n Rv VvgV, a.e. in (0,T) (2.35) 

99 ( 0 ) = cpo and cr(0) = cr 0 in ft (2.36) 

where R is defined by (I2.15p . accordingly. 


The above result generalizes the analogous [5J Thm. 2.6] as far as the assumptions on 
the initial data are concerned (and also since just the case a = f3 was considered there). 
Moreover, in [5j Thm. 2.6], even uniqueness for the solution to the limit problem was 
proved. However, also for this point, stronger conditions on the initial data are assumed 
in order that the solution to the limit problem is rather smooth. Here, we can consider 
the natural regularity requirements, i.e., 

h G T 2 (0, T; V) 

<p G H\ 0, T; V*) n L°°(0, T; V ) D L 2 (0, T; W) 
a G H\ 0, T■ V*) fl L 2 (0, T; V) C C°([0, T}; H). 


(2.37) 

(2.38) 

(2.39) 
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For uniqueness in this framework, we can quote the even more general result [T3J Thm. 2]. 
However, uniqueness also follows from the error estimate we present at once (see the 
forthcoming Remark H.ll for details). In order to state our last result we need to reinforce 
the assumptions we made on the potential F; namely, we assume 

D(B) = M and F = B + n is a C 2 function on K (2.40) 

\F(r)\ <C 0 (|r| 6 + l), \F'{r)\ < C^r] 5 + 1), and \F"(r)\ < C 2 (|r| 4 + 1). (2.41) 

Although the third condition in (12.41)1 implies the other two, we have written all of them 
for convenience. We also remark that the classical potential (12.8)1 fulfils such assumptions. 
Furthermore, we notice that (j2.41 jl is slightly more general than the analogous assumption 
made in j5, Thm. 2.6]. Finally, we can observe that the exponents in ()2.41jl are related to 
the dimension of R and the related Sobolev embeddings. Here is our last result. 

Theorem 2.7. Assume (12.31) - (12.7(1 and (12.40(1 (12.41)1 on the structure and (12.11(1 on the 
initial data. Then, with the notation of Theorem \2.(A the estimate 

\Wa,P — ^IU^( 0 ,T;V*)nL 2 ( 0 ,T;V) + || ha,/3 ~ 7t|U 2 (0,T;V*) 

+ || or o:,/3 — a \\L°°(o,T-,v*)nL 2 (o,T-,H) < C {a 1 ^ + fd 1 ^ 2 ) (2.42) 

holds true with a constant C that depends only on R, T, the structure of the system, and 
the norms of the initial data related to assumptions (12.1111 . but not on a nor on /d. 

The rest of the section is devoted to list some facts. We make repeated use of the 
notation 

Q t := fl x (0, t) for t 6 [0, T] (2.43) 

and of well-known inequalities, namely, of the elementary Young inequality 

ab < 5a 2 + b 2 for every a, b > 0 and 5 > 0 (2.44) 

as well as of Holder’s inequality and its consequences. Moreover, as fl is bounded and 
smooth, we can owe to the Poincare and Sobolev type inequalities, namely, 

IMIe < + | Jn v \) for every v G V (2.45) 

V C L g (Q) and ||u|| 9 < C||u||y for every v G V and 1 < q < 6 (2.46) 

L 9 (R) C V* and ||w||* < C||n|| g for every v & L q (Q) and q > 6/5 . (2.47) 

In (12.45(1 — ()2.47p . C only depends on R. Finally, we recall the interpolation inequality 

||n||^ < ||n||y ||u||* for every v G V, (2.48) 

which trivially follows from the identity ||n||^ = (v,v) for every v G V. 


3 Proofs of Theorems 12.21 and 12.3 


We start proving Theorem 12.31 in the following form: (I2.22p ((2.23(1 hold for every a and fd 
and every solution to problem (12.15(1 — (12.19)1 satisfying the regularity specified by (12.1211 - 
(12.14(1 . We do not know anything about well-posedness yet, indeed. 
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First a priori estimate. We test (12.161) and (12.181) by /i and a, respectively, and 
integrate over (0,t), where t G (0, T) is arbitrary. At the same time, we multiply (I2.17|) 
by — d t ip and integrate over Q t . Then, we add the resulting equalities to each other, 
obtaining 


? [ Ht)\ 2 + [ d t pp + [ |Vp| 2 

z Jn JQt JQt 

f ^dtp + p I \d t p\ 2 + \ [ \Vip{t)\ 2 + f F(if(t)) 

J Qt J Qt z Jn Jn 

+ A [ l cr ( i )| 2 + [ |Vcr| 2 + f R(a - p) 
z Jn J Qt JQ t 

= f / l/^°r + | / |V^o| 2 + [ F(<p o) + i [ |ao| 2 . 

Clearly, two terms cancel out. Moreover, F is nonnegative by assumptions (I2.5I) - (I2.6|) . 
Finally, we have R(a — p) = IS) 2 and |i?| < |S) sup ^Jp a.e. in Q with the notation (I2.2ip . 
Therefore, with the help of (12.41) we immediately deduce 

al ^ 2 |l/ i IU oo (0,T;H) + || ^ P\\l 2 {0,T;H) 

+ (3 1/2 \\d t <p\\L2(o,T-,H) + || V9?|| L oo( 0T;H ) + ||R(^)|li // i( 0i T;Li(n)) 

+ || ^|| L°°(0,T \H)C\L 2 (0,T\V) + 11 $ IU 2 (0,T;^) + II A|| L 2 (0,T;H) 

< C (<A 2 ll^olla + l|V^ 0 |U + ||F(«,)|l! / , 2 (n) + |k„|| H ) (3.1) 

for some constant C that depends only on p. Thus, in order to prove (I2.22|) . we have to 
complete the full norm of p and estimate the terms that are missing in (13.11) . 

Second a priori estimate. We estimate the mean value of p by testing (I2.20j) by 
v — 1. We obtain, for every t e [0, T], 

/ (ap(t) + p(t) + a(t)) = / (a/x 0 + p 0 + a 0 ) < |Q| l/2 \\ap 0 + <£o + &o\\h 
Jn Jn 

and deduce that (since a < 1) 


< pit ) < Cq^^WpoWh + ll^olln + ||cro||+ al ^\\^(t)\\H + ||°‘(^)||r/), (3.2) 

n 

where Cq depends only on fh 


Third a priori estimate. We test (12.18jl . written at the time t, with v(t), where v is 
arbitrary in L 2 (0, T; V). Then we integrate over (0, T ) with respect to t and obtain 


(d t a(t),v(t)) dt 


< (||Vcr|| L 2 ( 0) T;H) + \\R\\l^(0,T-,H))\\v\\l^(0,T-V) • 


This means that 


||^t cr IU a (0,T;V*) < ||Vcr|| L 2 (0,T;H) + || R || L 2 (0,T;H) ■ 


(3.3) 






















Colli 


Gilardi 


Rocca 


Sprekels 


10 


Fourth a priori estimate. Similarly, we test (12.201) . written at the time t, by v(t), 
where v is arbitrary in L 2 (0, T\V). We obtain 


{d t {ap + (p)(t),v(t)) dt 


< 


(d t a(t),v(t)) dt + ||V(/x + c)\\l 2 (o,t-,h)\\v\\l 2 (o,t-,v ), 


whence immediately 


\\d t {an + p)\\l 2 {o,t-,v*) < ll ^ cr l | L 2 ( o , T ; y *) + IIV/i||i2( 0T;H ) + || Va|| i 2(o i T;rr)- (3.4) 


First conclusion. We combine (I3.1IMI3.4I) with the Poincare inequality (I2.45P applied 
to and immediately deduce (j2.22j) with a constant C that depends only on p , R and T. 

Fifth a priori estimate and conclusion. By estimate (12.221) and the Lipschitz con¬ 
tinuity of 7 T, we deduce that 

ll 7r ( ( / 9 )IU 2 (o,T;jr) < C (a 1 / 2 ||/io||ir + ll^ollv + ll-f(^o)II za^o) + ll^olln + l) ; (3.5) 

with the same C , without loss of generality, provided that we allow C to depend on ir as 
well. Now, we write (I2.17|) in the form 

- A<p + f = / := -(3d t <p - vr((p) + p 

and observe that (I2.22p . (13.51) and (3 < 1 imply 

II/I|l 2 (o,t ; h) < C (oi 1 ^\\pq\\h + ll^ollv + ||^ r (v :, o)llzA^) + ll^olln + 1 + ||hlU 2 (o 

with the same C once more, without loss of generality. If M denotes the right-hand side 
of this inequality, a standard argument (formally multiply by —A<p) shows that both 
A<p and £ are bounded in L 2 {Q) by a multiple of M. Therefore, the same holds for 
IMU 2 (o,T;W) by elliptic regularity. Finally, the full norm ||hlU 2 (o,T ; v) is is equivalent to the 
sum of || V/i||L2(o,T;_ff) and ||A^-lli, 2 (o,cr ; Af)- Thus, (12.231) follows and the proof of Theorem 12.31 
is complete. □ 

Proof of Theorem 12.21 As far as uniqueness is concerned, we can refer to the proof 
of the uniqueness part of [5] Thm. 2.2] since it holds under the present assumptions. In 
order to prove the existence of a solution, we approximate the data /i 0 and cr 0 by functions 
Pq )£ and <to,£ satisfying 

ho.ej &o,e £ V for e > 0, po, £ —t ho and <to,£ —t cr 0 in H as e \ 0. 

Then, for every £ > 0, the condition (I2.24p holds for the approximating data so that the 
assumptions of [5J. Thm. 2.2] are fulfilled. Thus, the problem (12.15p (12.19p has a unique 
solution (/i £ , <p e , cx £ , £ e ) with R e defined by (12.15P accordingly. Moreover, such a solution 
must satisfy (I2.22MI2.23I) due to the above proof. As a and (3 are fixed, such estimates 
provide uniform boundedness with respect to e even for p £ in L°°(0, T; H) and d t p £ in 
L 2 (0,T; H). Therefore, (I2.23p implies that /i e , p £ and are bounded in L 2 (0,T;F), 
R 1 ( 0, T; H) flL 2 (0, T; W) and L 2 (0, T; H), respectively. Finally, the estimate for the time 
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derivative of ap £ + (p £ derived from ( 12 . 221 ) and the estimate for dtp £ mentioned before 
imply that d t p £ is bounded in L 2 (0,T; R*). Hence, we have 

/i £ -A n weakly star in R 1 (0, T ; V*) D L 2 (0, T; R) 

<p £ -A p weakly in if ^O, T; if) D L 2 (0, T; W) 

<r e —^ cr weakly star in if 1 (0, T; R*) D ii 2 (0, T ; R) 

£ e —>■ £ and R £ ^ R weakly in L 2 (0, T; if) 

as £ \ 0, at least for a subsequence. This implies, in particular, that the initial conditions 
for (p,, ip, a) are satisfied. Moreover, the above convergence for ip £ and the Aubin-Lions 
lemma (see, e.g., m Thm. 5.1, p. 58]) imply that 

fi £ —* /i, ip £ —* 99 , cr e —* cr strongly in L 2 ( 0 , T; if). 

Then, 7 r(</? e ) and p(yj e ) converge to 7 r(<^) and p(yj), respectively, strongly in L 2 (0,T;ff). 
Therefore, we can identify the limits of the nonlinear terms £ e and R e . For the former, we 
can apply, e.g., (T, Cor. 2.4, p. 41] and conclude that £ G B(ip) a.e. in Q. For the latter 
we note that R £ converges to p(ip)(a — /i) strongly in T 1 (Q), whence (j2.15[) follows. At 
this point, we can write the integrated-in-time version of problem (I2.16|) - (I2.18|) for the 
approximating solution with time dependent test functions and take the limit as e tends 
to zero. We obtain the analogous systems for (p, tp, a, £), and this implies (l2.16D - (l2.18j) 
for such a quadruplet. This completes the proof of Theorem 12.21 □ 


4 Asymptotics 

This section is devoted to the proof of Theorems 12.61 and 12.71 In order to simplify the no¬ 
tation, we follow a general rule in performing our a priori estimates. The small-case italic 
c without any subscript stands for different constants, that may only depend on R, T, the 
shape of the nonlinearities and the norms of the initial data related to assumption (12.111) . 
A notation like cs signals a constant that depends also on the parameter 5. We point out 
that c and c$ do not depend on a and /3 and that their meaning might change from line 
to line and even in the same chain of inequalities. On the contrary, those constants we 
need to refer to are always denoted by different symbols, e.g., by a capital letter. 

Proof of Theorem 12.61 We follow the argument done for [5] Thm. 2.6] rather closely, 
but we have to modify the types of convergence since our assumptions are different and 
more general. We start from (I2.22p - (I2.23|) . written for the solution <Pa,p, £a„a), 

and improve the latter by estimating the norm of on its right-hand side. However, 
we omit the subscripts a and (3 for a while. Thanks to (12.17)) and (12.261) - we have that 
l£l < C(B( V ) + 1) a.e. in Q. Then, by integrating over fl we obtain 

[ |£(t)| < C [ (B(<p(t)) + l) for a.a. t G (0, T). (4.1) 

Jn Jn 


At this point, we can estimate the mean value of /.! on account of (I2.6[) and (12.71) . Indeed, 
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by just integrating (12.171) over R, we deduce that 


A»(f) 


(Pd t <p + £ + n(ifi))(t) 


<flft»>Mlll+c(l|SWt))||i + Wt)||. + l) 

< c^\\a^p{t)\\ K + c(||F( V («))||, + IkWII* + 1) (4.2) 

for a.a. t G (0, T), beacause of the Lipschitz continuity of n and the nonnegativity of it. 
Then, (j2.22j) implies that the function t ha \f n p(t)dt\ is bounded in L 2 (0,T). By com¬ 
bining this with (I2.22p and the Poincare inequality (I2.45p . we derive that // is bounded 
in L 2 (0, T;V). Hence, recalling estimates (j2.22p (l2.23p it turns out that the convergences 
(I2.28p ~ (l2.32p and a convergence for R a p hold, at least for a subsequence. For the reader’s 
convenience, we write this conclusion explicitly, as well as the consequences we are inter¬ 
ested in. These are obtained by means of strong compactness results (see, e.g., [2H Sect. 8, 
Cor. 4]), the Sobolev inequality (I2.46P and the Lipschitz continuity of n and p. We have 

p a ,p —> /J weakly in L 2 { 0, T;V) fl L 2 (0, T; L 6 (f2)) 

p a ^ P weakly star in L°°(0, T;V) fl L 2 (0, T; W) 

v a ,p -> a weakly in H\ 0, T; V*) D L 2 (0, T; V) D L 2 (0, T; L 6 (R)) 

£,a,p —>■ i an d R e —^ R weakly in L 2 (0, T; H) 

ctRa,p —> 0 strongly in L°°(0, T; H) fl L 2 (0, T; V) D L 2 (0, T; L 6 (f2)) 

fidtPa,p —>■ 0 strongly in L 2 (0, T ; H) 

d t {ap a ,p + <Pa,p) -t dtp weakly in L 2 (0, T; V*) 

ap Q ,p + Pa,p P strongly in C°([0, T\; V*) fl L 2 (0, T; H) 

n{p a ,p) ->• n{p) and p(p a ,p) -A p(p) strongly in L 2 (0, T; H). 

Hence, we infer that p and a satisfy the initial conditions (12.361) . Moreover, we deduce 
that £ G B(p) (apply, e.g., p, Prop. 2.5, p. 27]) and that R a i( g also converges to p(p)(a—p) 
weakly in L 1 (0, T; L p (f2)) for some p G (1, 2): consequently, we have R = p(p)(a — p). 

Finally, we take the limit in the integrated-in-time version of problem (j2.16p ~ (l2.18p 
for {pa,piPa,pi a a,p,^a,p) with time-dependent test functions. We obtain the analogue 
for the system (12.33j) — (j2.35p . Finally, as the solution of the limit problem is unique by 
Theorem 14.11 the convergences we have obtained for a subsequence hold for the whole 
family. This completes the proof of Theorem 12.61 □ 


Proof of Theorem 12.71 As we use some ideas of |T3], it is convenient to rewrite 
the equations (12.16p and (12.181) as abstract equations in the framework of the Hilbert 
triplet (V, H, V*) related to an invertible operator. To this end, we introduce the Riesz 
isomorphism A : V —» V* associated to the standard scalar product of V, that is 

(Au, v) (u, v)v — / (Vu-Vu + uu) for u,vEV. (4.3) 

J n 

We notice that Au = —Au + u if u G W and that the restriction of A to W is an 
isomorphism from W onto H. We also remark that 

(Au,A _1 u*) = (v*,u) for every u G V and v* G V* 

{u*, A^v*) = ( u*,v *)* for every u*,v* G V*, 


(4.4) 

(4.5) 
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where (•, ■ )* is the dual scalar product in V* associated with the standard one in V, and 
recall that (v*,u) = f n v*u if v* G H. As a consequence of (14.5[) . we have 


d_ 

dt 



2 (d t v*,A V) for every v* G R 1 (0, T; V*). 


(4.6) 


In view of the regularity conditions (12.12p (12.14j) and (I2.37j) (j2.39j) . we rewrite (12.161) — 
(I2.18|) and (I2.33|) (I2.35j) for the solution (p a ,pi p a ,pi ®a,p) to (j2.15j) (12.19p and the one of 
the limit problem, respectively. If we term the latter we have 


& dtH'Oi “1“ dt(ficx.,f3 H - Afl^^ Ra,(3 “1“ 

(4.7) 

Mq,/5 “I - Ap a p F ( ^PQi^fS ) 3 

(4.8) 

dt&a,/3 H - A(J a,/3 Ra,/3 0"a,/3 

(4.9) 

dtp Ap = R fi 

(4.10) 

~jl = Alp + F'(jp) — Tp 

(4.11) 

d{o + Aa = —R + a, 

(4.12) 


where R a> g and R are defined by (12.15)) according to the equations we are considering. 
All these equations are meant in V* a.e. in (0, T). However, (14. 8 p and (14.lip also hold 
a.e. in Q. Moreover, the solutions have to satisfy the initial conditions (12.19p and fl2.3bH . 
respectively. Now, we take the differences between (14.71) (14.91) and (I4.10|) - (l4.12p at time 
s G (0, T ) and test them by 

A~ 1 (an at g + -(an at p + p)(s), and A _1 cr(s), 

respectively, where we have set for convenience 


H . 3 di P ■ Pa, ft Pi @ @ a, ft &i and R . Ra,g R- 

Next, we sum up and integrate over (0,f) with respect to s, for an arbitrary t G (0,T). 
We obtain (by omitting the evaluation at s inside integrals, for brevity) 


(d t (apt at g + </?), A 1 (a/j, at g + (p))ds+ / {Aft, A 1 (apt at g + </?)) ds 

J o 

- / (n,an a £ + (p)ds+ / {Pd t (p a> p,afjL at p + (p)ds+ / (A<p, afi a> g + <p) ds 

Jo Jo Jo 

+ / (. F\ip aj g) - F'(Tp), an a ,g + p) ds 


+ / (d t a,A 1 cr)ds+ / (Aa,A l a)ds 


= / (R,A (afj. a ,p + p)) ds + / {fi,A (apL^g + p)) ds + / {ip, api a ,p + </?) ds 


(R,A V) ds + / (a, A 1 a)ds. 


For the reader’s convenience, we just rearrange and use the decomposition F' = B + n. 
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We have 


/ (d t (a/i a ^ + (p),A 1 (a/i Q , i s + p)) ds 
Jo 

+ / (A/i, A~ 1 (a/j , a + p)) ds - / </x, a/jiafi + <yc) ds 
Jo Jo 

+ / (Ap,p)ds + / (B(p at p) - B(p),p) ds 
Jo Jo 

+ j (<9 t cx, dl _1 cr) ds + / (/ler,/l _1 cr) ds 
do do 

= / (-R, dl“ 1 (a/x ai/3 + p - a)) ds - / (/3d t (f a ^, a/i a ^ + p) ds - / (Ap, a/i a ^} ds 

Jo Jo Jo 

+ / (yii,dl _1 (a/x Q>/ 3 + p)) ds - / (7r(p aj/3 ) - n(p),p) ds 
Jo Jo 

- / (F'(p aj p) - F'(p),at/J. at p)ds+ / (<p, a/i a p + 99) ds + / (cr,/TV) ds. 
do do do 

At this point, we account for (14.3j) - (l4.61) and observe that the second and third terms on 
the left-hand side cancel out. Finally, owing to the initial conditions (ot/JL a ,p + p)( 0) = a/i 0 
and cr(0) = 0, we deduce 

^\\(<xfM Q ,f> + p)(t)\\l+ ^ \\p\\yds+ j (B(p atP )-B(p))p 

^ Qt 


^\wm: 


(7 


'Qt 


\ ||«/Ro||* + J (-R, a/ia t/ 3 + p - a)*ds - f {f3d t p a ,0,a/i a ^ + p) ds 
2 Jo Jo 

- / (p^a/i^p) ds + / (/x, ot/i a fi + p)*ds — / (ir(p a ^) - vr(^)) p 
Jo Jo JQt 


' Qt 


{F'(p a ,p) ~ F\p )) OL/lafi + / p(ap a p + p) + 

J Qt 


cr 


I? ds. 


(4.13) 


All of the terms on the left-hand side are nonnegative, the third one by monotonicity. Now, 
we treat each integral on the right-hand side, separately. In the sequel, 5 is a positive 
parameter whose value will be chosen at the end of the procedure. We first observe that 
(12.22)) holds for the solution (fi a ,Pi Va,p, and that Theorem 12.61 improves (12.231) for 
such a solution. Indeed, the restricted setting of regular potentials satisfying (12.41 jl led 
to (14.21) . So, as we have seen in the previous proof, (12.22)) and (12.23)) imply 


\\da,p\\L 2 (0,T-y) + \\Poi,p\\l 2 (0,T-,W) < c. (4-14) 

Now, we prepare estimates for ||/i||* and ||i?||* a.e. in (0, T). Again for simplicity, in 
performing them, we omit writing the evaluation point. From the mean value theorem 
and the third assumption in (12.411) we easily derive that 

\F\p a ,p) - F\p )| < c\p\(\p a p\ A + \p\ 4 + 1) a.e. in Q. 
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Therefore, by the Holder and Sobolev inequalities, we infer that 
II F '( ip at p ) - F '( jp )\\ 6 / 5 < C 11 ^ 116 ( 11 ^ 0 ,^ 113/2 + ||^ 4 || 3 /2 + l) 

= C MUdl^lli + Irt + 1 ) < cM\v(\M\v + m\v + 1 ) < c||<p||v, (4.15) 

the last inequality following from estimate (12.22 1) for ip a ^ and the regularity d2.38[i of Tp. 
Taking the difference between (14. 8p and (14.111) and using the dual Sobolev inequal¬ 
ity (j2.47p . we deduce that 


IImII* = \\Pdt<Pa,l3 + A(p + F'iVafi) - F\Tp) ~ <f\\* 

< P\\dt<p a ,l3\\* + |M|v + c\\F\(p aj p) - F\Tp) || 6/5 + IMI* 

< c{p\\dtip a ,fi\\* + |M|y) < c/3 1 / 2 + c \\<p\\ v , (4.16) 

where the last inequality follows from (12.221) . In order to estimate ||i?||*, we first observe 
that the boundedness and the Lipschitz continuity of p and the Sobolev inequality (applied 
to Vp and the test function v E V) imply that, for every v E V, 

\\p(jp)v\\ v < \\p{Jp)v\\ H + || Vp(^) v|| H + \\p(jp)Vv\\ H 
< c\\v\\h + c||V^|| 4 IMU + c\\\7v\\ H < c(||^||w + l)|M|v • 


Hence, we have for every v E V the estimate 


Rv= (p((p aiP )(a a} p - p ai p) - p(ip)(a - fj))v 

J u J u 

< / b(<A*,/0 -p(jp)\ Wa,p-p a ,M M + / p(jp)(a-p)v 

Jn Jn 

< c||<p|| 3 \\a a> p - Pa^UMU + lb - p\\* \\p<&) v\\v 

< c\\ip\\v \\o a ,p - ^llvlblb + c|b - bl* (WWw + l)\\v\\v 


Therefore, we can estimate ||i?||* a.e. in (0, T), also owing to (I4.16p . in this way: 


||R||* < c|Mb \Wa,P - Pa,/3\\v + c\\<T ~ p ||* (|Mlw + 1) 

< C|b|b || (T at p - Pa,MW + c(|b||* +/3 1/2 + |M|y) (||^||lV + 1) 

< C^aA |Mlv + INI*) +C/3 1/2 -0, 


where ba,/ 3 > b : (0, T) > M. are defined by 

'Ma,?'■= \Wa,p - Pa,p\\v + W\\w + 1 and b := ||^lk + 1, a.e. in (0 ,T). 
Coming back to the right-hand side of (14.131) . we can treat the first term as follows: 


/ (R, ap^ + (p - a)* ds < / ||-R||* \\ap a> p + ip « cr||* ds 
Jo Jo 

< / (cba,/l(|b||y + ||(t||*) + C/3 1/2 b) (ll«4ta,/3 +/2||* + ||a||*) ds 

Jo 

<5 \\tp\\v ds + (3 [ |b| 2 ds + c s [ 'ipl, 0 {\\aPa,p + <p\\l + \W\\l)ds 

Jo Jo Jo 


(4.17) 
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We observe at once that the regularity 02.381) for Tp and estimates 02.22|) for <j Q ^ and 01.11 1) 
for imply that 0 G L 2 (0,T) and that -0 aj p is bounded in L 2 (0,T), so that 0 2 i/5 is 
bounded in L 1 (0, T). This will allow us to apply the Gronwall lemma. Now, we estimate 
the next term on the right-hand side of 04.131) . Using (12.22(1 . we see that 


- / {/3d t ip a ^,ap a}0 + (p)ds 

Jo 

— a II/G^IIa 2 (o,T ; I?) + ^ / IMItf ds + Cs j3 \\dt<fia,/3\\L 2 (0,T; H ) 


< ca + 8 


V 


ds + cs (3 ■ 


Next, we have 


(4.18) 


(4.19) 


- / (ip,ap atf} ) v ds < 8 / \\ip\\y ds + c s a 2 
Jo Jo 

thanks to 04.14|) for p, ai0 , as well as, by 04.16|) . 

/ (/.i,a/j, aj p + ip)*ds <c (/? 1/2 + IMIv) \\oi(J>a,p + V?||* ds 
Jo Jo 

<8 [ \\v\\v + P + C8 [ \\aiJ. a ,p + <p\\Ids. (4.20) 

Jo Jo 

Moreover, by using the Lipschitz continuity of n, the interpolation inequality 02.481) and 
(14.14)1 for once more, we can write 


' Qt 


{n(<Pa,p) - ir(<p)) V < c \<p\ 2 < c / \\(p\\v 

J Qt Jo 


; dS 


< 8 


yds + cs 


: ds 


< 5 


yds + cs / ||a/z Q ,/j + v\\lds + c 5 a 2 . 


(4.21) 


The next term to deal with is the one involving F'. We use 04.1511 . the Holder, Sobolev 
and Young inequalities, and the estimate 04.141) for Thus, we have 


1 Qt 


(F'(<Pc,f)) - F'(<p)) ap aJ 3 < / II F\y a ,p) - F'(ip )|| 6/5 ||a/x ai ^|| 6 ds 


< c / \\(p\\ v a\\p at p\\ v ds < 5 


y ds + cs or 


(4.22) 


Finally, the last integral on the right-hand side of 04.13)1 does not need any treatment and 
the preceding term can be estimated in this way: 


(<£>, Otp a fi + ip) ds < 8 


yds + cs / ||a/x a>/3 + ^||*ds. 


(4.23) 


At this point, we combine (14.13)) and the list 04.171) 04.23)) of estimates we have obtained. 
Then, we choose 8 small enough, recall that 0 G L 2 (0,T) and that 0 2 i/3 is bounded 
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in L 1 (0,T), and apply the Gronwall lemma in the form |2J Lemma A.4, p. 156]. We 
obtain 

2 ll( a /W£ + tOWII* + J IMIv ds + 2 IkWII* + J M 2 < c ( a + / 5 ) 

for every t G [0, T\. As ||ct:/i a , i/ 3 (t)|| 2 < ca for every t G [0, T] by (12.221) . the above inequality 
implies 

IMlL^( 0 ,T;V*)nL 2 ( 0 ,T;V) + II^llL°°( 0 ,T;V*)nL 2 ( 0 ,T;H) < C + Z? 1 ^ 2 ). (4.24) 

Now, we take the differences of equations (14.81) and (14.111) and estimate the F 2 (0,T; V*) 
norm of it. With the help of (I2.22|) and (I4.15|) it is straightforward to infer that 

IHIl 2 ( 0 ,T;V*) < c/3 \\d t ip a fi\\ L 2(Q'T-,H) + C || <~P || L 2 ( 0 ,T; V) 

+ c II F'(ip Qi p) - F'(jp) ||x,2( 0)T;L 6/5( r2 )) + |MU 2 (o,T;V*) 

< c/3 1//2 + c ||</3||L2( 0i T;y)- (4-25) 

Hence, in view of (I4.24p and (14.251) we finally obtain the estimate (I2.42p . where one has 
to read Tp, JI and a in place of p>, p and a, respectively, due to the change of notations 
within this proof. □ 

Remark 4.1. By going through the above proof, one sees that uniqueness for the limit 
problem (I2.33p - (l2.36p has been never used, that is, the following formulation of Theo¬ 
rem EH has been proved: the error estimate (I2.42p holds for every solution (/u, ip, a) of 
the limit problem satisfying the regularity requirements (12.371) (12,39j) . This implies the 
uniqueness for such a solution. Indeed, if (/ij, <pi, cp), i = 1,2, are two solutions of the 
limit problem, by writing (12.421) for both of them and using uniqueness for the solution 
{p a p, H a ,p, cr a ,p)i one immediately derives 

IIVT — 7^2 \\ L°°(0,T-,V*)nL 2 (0,T,V) + ||hl — h21| L 2 (0,T;V*) 

+ ||(Ti — 0"21|L°°(0,T;V*)nL 2 (0,T \H) < C (d? 1 ^ 2 + 

for every a, (3 E (0,1), whence p\ = tp 2 , Hi = /i 2 and 0 \ = ct 2 . Then, by comparison in 
(I2.34p . it follows that £1 = £ 2 , as well. 
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